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Summary: Let a > 0, 8 > 1, b > 0, < p < 1. Let A be a continuous and positive 
function in H^(M. + ). Using the technique of moving domains (see [14]), and classical 
direct stochastic calculus, we construct for positive initial conditions a pair of continuous 
and positive semimartingales (R, yR) with 



> ' dRt = o\[RtdW t + —(8- bR t )dt + (2p - l)d£° t (R - A 2 ) 

0> ! 4 

and 



dVR t = ^-dW t + ^ ( 5 —=^ - bVRt] dt + (2p - l)df t (VM - A) 
2 8 V VRt J 

+ I -^£° t + (VR), 



where the symmetric local times £°(R — X 2 ), £°(y/R — A), of the respective semimartingales 
R — A 2 , \/R — A are related through the formula 



2VRd£°(VR - A) = d£°(R - A 2 ). 

Well-known special cases are the (squared) Bessel processes (choose a = 2, b = 0, and 
A 2 = 0, or equivalently p — |), and the Cox- Ingersoll- Ross process (i.e. R, with A 2 = 0, or 
equivalently p = |). The case < 5 < 1 can also be handled, but is different. If \p\ > 1, 
then there is no solution. 
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1 Introduction 



Consider a continuous, and positive function A : R + — > R + , A G For parameters 

cr > 0,5, b > 0, p G]0, 1[, and positive initial conditions, we construct the Cox-Ingersoll- 
Ross process with skew reflection on A 2 (here X 2 (t) = X(t) ■ X(t)), i.e. a weak solution 
to 

2 

dR t = ay/\R^\dW t + ^-(<J - bR t )dt + (2p - l)df t (R - A 2 ), (1) 

where £^(R — X 2 ) is the symmetric local time at zero of the continuous semimartingale R — 
A 2 . The equation is in particular complicated because of the singular diffusion coefficient 
in combination with a parabolic local time. For example, the first naive idea to put 
S t = Rt — A 2 (t) in order to simplify (1) doesn't improve at all the situation. One ends up 
with 

dS t = ay/\S t + \*(t)\dW t + ^(5-bS t - b\ 2 (t) + ^X(t)X'(t))dt + (2p - l)df t (S), 

4 a z 

i.e. an equation with degenerate time dependent diffusion coefficient and still remaining 
reflection term. For equations of this type general existence and uniqueness results are 
completely unknown. Playing around with equation (1) one quickly observes that the sole 
employment of direct stochastic calculus doesn't seem to be helpful for the construction of 
such highly singular processes. In fact, constructing a solution to (1) means to construct 
a solution to the Cauchy problem related to the generator of the CIR process with skew 
boundary conditions on A 2 . This is in view of the coefficients a complicated task. Our 
approach here is to provide a general method for the construction of these type of diffusions 
applying a recently developed technique of moving domains (see [14]) in combination with 
classical direct stochastic calculus. Let us explain this method in detail. We keep 8 > 1 
in order to simplify (for the case 5 G (0, 1) see Remark 2.9). First one can readily check 
that (1) is the square of a solution Y, Y > to 

dY t = Uw t + y - bY^j dt + (2p - l)I {A>0} ^ t (F - A) + Ii|i>^+(y), (2) 

where the upper (right-continuous) local time. In fact, just apply Ito's formula to see 
that R t := Y t 2 , R > 0, solves 

2 

dR t = a^RtdWt + bR t )dt + (2p - l)2y^d£° t (VR - A), 

and then notice that 

2y/R t d£° t (v r R - A) = d£° t (R - A 2 ). (3) 

The relation (3) can be shown by probabilistic means using a product formula for local 
times (see [18], and [8]). In fact, we have R - A 2 = {y/R + A)(v^R - A), but see also 
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Remark 2.8(i) for an independent analytic proof. So far we explained why it is enough to 
construct (2) in order to get (1). The advantage of (2) in comparison with (1) is that it 
is much better suited to be obtained by a change of measure since there is no coefficient 
in the martingale part. 

The next step is to decompose A into A = (5 + 7, where (3, 7 G H^(R + ) are continu- 
ous representatives, with f3 decreasing, and 7 increasing. Consider the following moving 
domain 

E := {(t,x) G R + x R\x > -7(f)}. 

Since —j(t) decreases, E t := {x G rr) G E} = [—7(^,00) is increasing in t. As 
reference measure on E consider 

m{dxdt) := p(t,x)dxdt, 

where 

p(t,x) := ((1 - p)I { ^ mm) {x) + P I [moo) (x)) \x + ^(t)\ s - l e- h -^- 
is assumed to be increasing in t, that is 

p(s, x) < p(t, x) V < s <t, x G E s . 

Due to the monotonicity assumptions on p (and E) we are in the theoretical framework of 
[14]. Note that the monotonicity assumptions are in fact restricting as these do not allow 
all imaginable choices of parameters. We will comment on this below in the introduction. 
Applying the theory of time-dependent Dirichlet forms on monotonely moving domains 
of [14] we deduce that the process associated to the time-dependent Dirichlet form (10) 
solves 

X t = X Q + ^B t + ^ f 5 ~ l -bX s ds 
2 8 Jo X s + 7(s) 

+ (2p - 1) jy W s)>-,(s)}dt s (X + l -^^t + {X + 7 ), (4) 

(that is (28)). The construction is carried out in section 2 with advanced and heavy 
machinery, i.e. the Theory of Generalized Dirichlet Forms. We couldn't see how else to 
do it. In fact there is a lot of work done with direct stochastic calculus in the case of 
"nice" coefficients and non-parabolic local time (see e.g. [5] as showcase), or also in the 
Brownian motion case (no coefficients) with parabolic local time (see e.g. [2], or [16]). 
So, the type of equation (1) seems to be a real novelty. Now let us come back to the 
explanation of its construction. In section 2.1.1 we first derive the existence of a diffusion 
associated to the time-dependent Dirichlet form (10) (see Theorem 2.2). In order to obtain 
finally (4) we first identify (22) as the process associated to the time-dependent Dirichlet 
form (10). Since in general Dirichlet form theory the drift is associated to a signed smooth 
measure via the Revuz correspondence (14), we can see Markovian local times, i.e. singular 
increasing processes that are associated to smooth measures and their potentials, in the 
equation (22) (cf. explanation at the beginning of section 2.1.3). That these Markovian 
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local times exist is shown in Proposition 2.3. In particular the signed smooth measure 
corresponding to the Markovian local times in (22) is given in (17) with F equal to 
the identity (see also (19), (20)). The Markovian local times are indeed semimartingale 
local times. This is deduced in section 2.1.3 by comparing Fukushima's decomposition for 
generalized Dirichlet forms (16) with Tanaka's formula (23) (see explanation right after 
Remark 2.4). The final formulas in order to obtain (4) are then given in (26), and (27). 
Once a solution to (4) is constructed we just define 

1 r* 

W t :=B t + — 8 7 '(s) + <x 2 6 7 ( S )ds, 
4(T Jo 

and then under a change of measure with density 

e -h Jo 8l'(s)+^l n (s)dB s -^ 7 J* |8 7 'W+ CT 2 b 7 ( S )| 2 d S 

we see that Y t :— X t + 7 (i) is a solution to (2), hence R := Y 2 solves (1) and R is con- 
structed. 

Let us now comment on the monotonicity assumptions on p (and E). The time dependent 
Dirichlet form (10) consists of two parts. The first is a symmetric bilinear form and the 
second is a perturbation with the time derivative. The time derivative is generated by the 
semigroup (U t )t>o of (9). U t pushes the support of a function F G Cq(E) by t to the left. 
If we want to guarantee that (U t )t>o becomes a contraction semigroup on L 2 (E , pdxdt) 
then this is only possible if the support of UtF is smaller than the one of F (i.e. E t in- 
creases in t), and if p(s — t,x) is smaller than p(s,x) (p increases in t). The semigroups 
corresponding to the symmetric part of the time dependent Dirichlet form (10) and to 
the time derivative are nearly "orthogonal". Therefore in order to guarantee that the 
generator corresponding to the whole form (10) satisfies the positive maximum principle 
(cf. e.g. [4, chapter 4]) we need (U t )t>o to be a contraction. Various considerations lead to 
the belief that this might be difficult if not impossible to improve. In any case this is an 
intriguing question. For example it is not possible to add some killing with a constant c to 
make (e~ ct U t )t>o a contraction (and then to remove afterwards the killing). If p G (|, 1), 
so that 1 — p < p, then p(-,x) always increases, or if j3 — const, then p(-,x) increases 
for any p G (0, 1). Therefore our construction takes place for arbitrary A G Hj^(R + ), if 
2p — 1 > 0, or for increasing A, if 2p — 1 < (cf. below (6), and Remark 2.4, 2.8(h) for 
more general A). It is remarkable that pathwise uniqueness in [15] could be deduced for 
2p — 1 < if (A 2 )' > ^(5 — b\ 2 ). Thus for p < | and some strictly decreasing A but 
which is still bounded below by the mean-reverting level |. We also find remarkable that 
we were not able to construct a solution in the extreme cases p = and p = 1. 
In Remark 2.7(i) we show that if \p\ > 1, then there is no solution to (1). The construction 
of R, y/~R is in the sense of equivalence of additive functionals of Markov processes (see 
Remark 2.6). The case < 5 < 1 can also be handled, but is different, see Remark 2.9. 
In [15] it is shown that a solution to (1) always stays positive when started with positive 
initial condition. One can hence discard the absolute value under the square root in (1) as 
in the classical situation where A 2 = 0. Nota bene, that in case A 2 = 0, the local time £(R) 
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vanishes as a direct consequence of the occupation time formula. Let us also remark that 
the parallel work [15] is completely different from this work as well as from the subject as 
from the used techniques since there exclusively probabilistic arguments are employed. 
We really think that our approach is a novelty, and that equations (1), and (2) are truly 
worth to be studied. The parallel work [15] which deals with pathwise uniqueness of (1) 
shows that stochastic calculus is possible, even with such highly singular equations. It 
encourages to go further. For instance it could be challenging to further investigate er- 
godic behavior of (1) under the assumption lim^ 00 A 2 (t) = const. Or, is it possible to 
write down the distribution of R and/or to calculate explicitly the relevant correspond- 
ing quantities (mean, variance, etc.)? Finally we think that equation (1) can also have 
suitable interpretation in terms of multi-factor term structure models (i.e. models corre- 
sponding to a system of SDEs of basic univariate models) if we manage to provide enough 
analytic tractability and a sufficiently general framework to express no-arbitrage in the 
corresponding model. 



2 Construction of the skew reflected CIR process 
2.1 Construction of a solution to (4) 

2.1.1 The Generalized Dirichlet Form associated to a solution of (4) 

Throughout this article Ia will denote the indicator function of a set A. Let E := R + xR + , 
where R + := {x G R\ x > 0}. Let C^(R+) = {/ : R + -> R\ 3u G Cq(R) with uI K + = /}, 
and Cq(M) denotes the continuously differentiable functions with compact support in 
R. Let H 1 ' 2 (R + ) be the Sobolev space of order one in L 2 (M + ), that is the completion 
of Cq(1R + ) w.r.t. |0|ffi,2(R) = (J R+ |<9«0| 2 + \(fi\ 2 du)2. When considering an element of 
H 1 ' 2 (R + ), we always assume that is it continuous by choosing such a version. Later we 
will use the notions d u , du for the space variable (notation u = x) as well as for the 
time variable (notation u — t, or u — s). For the space-time variable we use y, e.g. 
y = (s,x), y = (t,x). Let H^(R + ) denote the space of all continuous : R + — > R such 
that (f)f E iJ 1 ' 2 (M+) for any / G C 1 (M+). 

If A G Hl^(R + ), then it has a uniquely determined continuous version w.r.t. the Lebesgue 
measure dt. We will always assume that A is continuous. Furthermore we assume that A is 
positive, i.e. A > 0. In particular A = /3 + 7, where f3, 7 G Hf ( £(R + ), and j3 is decreasing, 7 
is increasing. Indeed, since d t X G L 2 oc (IR + ) we may consider its positive part (<9 t A) + , and 
its negative part (<9 t A)~, and fix from now on 

/3(t) :=- [ (d t X)-(s)ds + X(0), 7 (t):= [ (d t X) + (s)ds. (5) 
Jo Jo 
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A = (3 + 7 > implies (3 > —7. Consider the following moving domain 

E := {(t,x) 6l + x R\x > ~7(t)}. 

Observe, that its t-section E t = {x G M.| (t, rc) E E} = [—7(^,00) is increasing in t since 
— 7(t) decreases in t. In particular A = U t >o{t} x E t . 

Let 5 > 1, b G M + (for the case <5 G (0, 1) see Remark 2.9). As reference measure on E we 
take 

m(dy) = m(dxdt) := p(t,x)dxdt, 

where 

p(t,x) := ((1 -p)I himm) (x) +pI m ,oo)(x)) |a; + 7(t)|^ 1 e-^ 
is assumed to be increasing in t, that is 

p(s, x) < pit, x) V < s <t, x G E s . (6) 

For instance, if p G (|, 1), so that 1 — p < p, then p(-, x) always increases, or if (3 = const, 
then p(-,x) increases for any p G (0, 1). 

Due to the monotonicity properties of p (and E) we are in the framework of [14]. More 
precisely, there is a time-dependent generalized Dirichlet form 8 with domain jFx VUVxjF 
on Ti := L 2 (E,m) which we determine right below. For q > 1 let 

C q (E) := {/ : E -> R\ 3u G C 9 (M 2 ) with ul s = /}, 

and Cq{R 2 ) denotes the g-times continuously differentiable functions with compact sup- 
port in R 2 . Let < o G R 

2 roo poo 

A(F,G):=— / d x F(s,x)d x G(s,x)p(s,x)dxds; F,GeC l (E), (7) 

o Jo J--y(s) 

with closure (^4, V) in H. The closability easily follows since p satisfies a Hamza type 
condition (see [14, Lemma 1.1]). Let A a (F, G) := A(F, F) + a(F,F), a > 0, where (•, •) 
is the inner product in Ti. For K C E compact, the capacity related to A is defined by 

Cap^AO = inf{A(A,F);F G C^E)}, (8) 

where Cq K (E) = {F G C^(E)\F(s, x) > l,V(s,x) G K}. For general A C A it is extended 
by inner regularity. Define 

U t F(s,x):= F(s + t,x); F E C l (E). (9) 

It then follows from results in [14] that (U t )t>o can be extended to a C -semigroup of 
contractions on TL which can be restricted to a Co-semigroup on V. For the corresponding 
generator (d t , D(d t ,H) on H it follows that 

d t :D(d t ,nnv)^V 
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is closable as operator from V to its dual V (see [10, I. Lemma 2.3.]). Let (d t , J 7 ) be the 
closure. T is a real Hilbert space with norm 



\F\t:= sJ\F\%+\d t F\l,. 

The adjoint semigroup (U t )t>o of (U t )t>o in H can be extended to a Co-semigroup on V'. 
The corresponding generator (A, -D(A, V')) is the dual operator of (d t , D(d t ,V)). T : = 
.D(A, V) fl V is a real Hilbert space with norm 



\F\ r -=W\%^\^F\V 

Let (•, •) be the dualization between V' and V. The time- dependent generalized Dirichlet 
form is now given through 



£(F,G) := 



A(F, G) - (d t F, G) for F e T, G e V 
A(F, G) - (AG, F) for Gef, F G V. 



Note that (•, •) when restricted to 7i x V coincides with the inner product (•,•) in 7i. 
In particular when F G Cq(F), G G V, then 

^2 poo poo 

£(F,G) = — / / d x F(s,x) d x G(s,x)p(s,x)dxds 

8 j y_ 7 ( s ) 

— / / d t F(s,x)G(s,x)p(s,x)dxds. (10) 

io V-7(s) 

For all corresponding objects to £ which might not rigorously be defined here we refer to 
[14] . We also point out that the monotonicity assumption on E t as well as on the density 
p in time is crucial for the construction of £. 

By [14, Lemma 1.6, Lemma 1.7] the resolvent (G a ) a>0 and the coresolvent (G a ) a>0 associ- 
ated with £ are sub-Markovian and Cq(E) C T dense. Let £ a (F, G) := £(F, G)+a(F, G) 
for a > 0. Then 

£a(G a F, G) = (F, G) H = £ a (F, G a G) F,G eV. 
Proposition 2.1 (G a ) a>0 is Markovian, i.e. GiI E = Ie m-a.e. 

Proof In order to prove the conservativity of £ it is enough to show that for one F G 
TL D L\E,m), F > m-a.e., there exists (W n ) n >i C < W n < I E , n > 1, W n | I E as 
n — > oo, such that 

lim £(W n ,G 1 F) = 0. 

Indeed, if this is the case then 

= lim £{W n ,G 1 F) = lim / (W n - G x W n )Fpdxds = I (I E - G^Fpdxds, 

' E J E 



n—*oo 
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and GiI E = Ie as desired. We now fix F as above, and determine below (W n ) n >i. 

Let g n e C 1 (1R + ), u n G Cq(R), n > 1, such that < g n ,u n < 1, \d t g n \oo, |<9 x m„|oo < L-n' 1 , 

\dxxUn\oo < L ■ n~ 2 , where L is some positive constant, and 



9n(s) 

and 



1 if s e [0, n\ 
if s e [2n, oo), 



u„(x) 



1 if [-7(2n)] < x < [A(0) + 1] + n 
if x > [A(0) + 1] + 2n, 



where [x] := sup{A; e 1\k < x}. Then VF n := g n u n ^-E € JF, n > 1, satisfies W n j as 
n — > oo, and since 9 x W n (s, /5(s)) = c^W^s, — 7(s)) = for all s, we easily find 

2 /*oo ;>oo / /A 1 \ 8 \ 

S(W n , G 1 F) = " y ^ y } [9 xx W n + y (x + ^ s)) - bx ) d * w n + —A W n) GiFpdxds, 
so that \S(W n ,GiF)\ is dominated by 

8 /> i[X(0)+i]+„ V'! 2 n([A(0) + l]+n + 7 (0)) n ) 



p2n poo 

+ / —GxFpdxds. 
Jo J 



'0 J-y(s) 



n 



Noting that GiFpdxds is a finite measure, we just apply Lebesgue's theorem, and the last 
sum is easily seen to converge to zero as n — > oo. This concludes the proof. 

□ 

Let us define the strict capacity corresponding to S. We fix $ e L 1 (i?,m), < $ < 1. 
Let A l)u be the 1-reduced function of kG^ A 1 := minftG^, 1) on [/, and let 

Cap 1 g 1$ (C/) = lim / A l)[/$dm if C/ C £ is open. 

If A C -E arbitrary then 

Ca Pl g iS (A) = inf{Cap 1; g i$ (C/)|C/ DA,[/ open}. 



We adjoin an extra point A to E and let £a := -EU{A} be the one point compactification 
of E. As usual any function defined on E is extended to E& putting /(A) = 0. Given an 
increasing sequence (F k ) keN of closed subsets of E, we define 

CoodFk}) = {/ : A -> R | (J F fc C A C -E, /|F fc u{A} «s continuous Wk}. 

k>i 
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A subset N C E is called strictly ^-exceptional if Cap 1( g ^(iV) = 0. An increasing se- 
quence (-F fc ) fceN of closed subsets of E is called a strict £-nest if Cap : q $(F%) { as 
k — > oo. A property of points in E holds strictly £-quasi-every where (s.S-q.e.) if the 
property holds outside some strictly ^-exceptional set. A function / defined up to some 
strictly ^-exceptional set iV C E is called strictly ^-quasi-continuous (s.£-q.c.) if there 
exists a strict £-nest (F k ) k€N , such that / G C OG ({F fe }). 

For a subset Ac£ A let a A := inf{t > | Y t G A} (resp. D A = inf{t > 0|F t G A}) be 
the /irsi hitting time (resp. /irs£ entry time) w.r.t. M. For a Borel measure u on E and a 
Borel set -B let P V (B) := f E P y (B)u(dy) and .E^ be the expectation w.r.t. P v . As usual we 
denote by E y the expectation w.r.t. P y . If U C E is open, then 

Cap 1>9l# (tf) = / S s [e-^]$(y)m(dy). (11) 

If £? C -B is an arbitrary Borel measurable set, then 

Cap lA# (S) = / ^[e- D -]$(y)m(dy). 

Both follows from [13, Lemma 0.8]. 

By strict quasi-regularity every element in T admits a strictly £-q.c. m-version (see [13, 
Proposition 0.9]). For a subset V cH denote by V all the s.£-q.c. m-versions of elements 
in T>. In particular denotes the set of all s.£-q.c. pdy- versions of 1-excessive elements in 
V which are dominated by elements of T . We have an analogy, namely [13, Theorem 0.16], 
to [11, Theorem 2.3]. That is: Let u G Pp. Then there exists a unique a-finite and positive 
measure on (E, 13(E)) charging no strictly ^-exceptional set, such that 

/ f dfiu= lim £i(f,aG a+ iu) V/' G Pp — V? . 
Je a ^°° 

Also in analogy to [11] we introduce the following class of measures 

5oo := {//a | u G Vq iH + and < °°} 

where G{Ht ■= {&ih \ h G 

For B G 13(E) the following is known from [13, Theorem 0.17]: B is strictly ^-exceptional 
if, and only if fi(B) = for all \i in S'oo- 

Since (£,T) is regular, i.e. Cq(E) n T is dense in Cq(E) w.r.t. the uniform norm as well 
as in J 7 , it follows that (£, T) is a (strictly) quasi- regular generalized Dirichlet form on E. 
On the other hand we can find a dense algebra of functions, namely Cq(E), in T . These 
two facts imply the existence of a Hunt process associated to £. Applying additionally 
Proposition 2.1, and [14, Theorem 1.9] we have: 

Theorem 2.2 There exists a Hunt process M = (Q, (J 7 t)t>0j (Yt)t>o, (Py)y=(s,x)eE A ) with 
state space E, and infinite life time, such that R a F(s, x) := J Q e~ at F {Y t (u)) P^ s>x )(dui)dt 
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is a S-q.c. m-version of G a F for any a > and any F G Tib- Moreover there exists a 
£ -exceptional set N C E such that 

P(s,x) {t | — ► Y t is continuous on [0, oo)) = 1 for every (s,x) G E\N. 
We will see that M = (Q, (.F t ) t >o, (Y t )t>o, (Py) y =(s,x)eE A ) in fact solves (4). 



2.1.2 Identification of the process Y associated to the Generalized Dirichlet 
form involving the Markovian local time 



We want to identify Y. Let us first recall some basic definitions and facts about additive 
functionals related to generalized Dirichlet forms. 

A family (A t ) t > of extended real valued functions on SI is called an additive functional 
(abbreviated AF) of M = (SI, {F t )t>o, (F t ) t > , (P y ) y eE A ) (w.r.t. Cap lA J, if: 

(i) A t (-) is ^-measurable for all t > 0. 

(ii) There exists a defining set A G J^oo and a strictly ^-exceptional set N C E, such 
that Pj,(A) = 1 for all y G E \ N, t (A) C A for all t > and for each uo G A, t h-> AM 
is right continuous on [0, oo) and has left limits on (0, CM), -AqM = 0, |A t M| < oo for 

t < CM, AM = AM for * ^ CM and A+ S M = AM + A(#M for s,t > o. 

An AF A is called a continuous additive functional (abbreviated CAF), if t i— > AM 
is continuous on [0, oo), a positive, continuous additive functional (abbreviated PCAF) if 
AM > and a finite AF, if | A t (uj) |< oo for alU > 0,u G A. Two AF's A,B are said 
to be equivalent (in notation A = B) if for each t > P y {A t = B t ) = 1 for strictly £-q.e. 
y £ E. The energy of an AF A of M is defined by 

e(A) = lim ^a 2 ^pdy 

whenever this limit exists in [0, oo]. We will set e(A) for the same expression but with lim 
instead of lim. 



-at \2 



Afdt 



(12) 



Let F be a strictly £-q.c. pdy-version of some element in Ti. The additive functional 

A^ := (F(Y t ) - F(Y )) t >o 

is independent of the choice of F (i.e. defines the same equivalence class of AF's for any 
strictly £-q.c. pdy- version F of F). The sub-Markovianity of (G a ) a>0 implies 

e(A^) = TmT (a(F — aG a F, F) H — ^ / (F 2 — aG a F 2 )pdy\ 
< Tim" a(F — aG a F, F) n . 
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Since T C V T (cf. e.g. proof of [12, Lemma 3.1]) it follows lim^oo aG a F = F weakly in 
V. Hence lim a ^ 00 a(F - aG a F,F) n = \im a ^ 00 £(F,aG a F) = S(F,F) whenever FGf. 
In particular 

e(A [F] ) < 2\F\ 2 r foranyFG.? 7 . (13) 

Define 

M = {M\M is a finite AF, E y [M?} < oo, E y [M t ] = 
for strictly S-q.e y G E and alH > 0}. 

M G M. is called a martingale additive functional (MAF). Furthermore define 

M = {M G M\ e(M) < oo}. 

o 

The elements of Ji4 are called MAF's of finite energy. 

Let A be a PCAF of M. Its Revuz measure ha (see [11, Theorem 3.1]) is defined by 



/ G(y)fi A (dy) = lim aE pdy [ e - at G(Y t )dA t 



for all G G B + . (14) 



The dual predictable projection (M) of the square bracket of M EM. is a PCAF of M. It 
then follows from (12), (14), that one half of the total mass of the Revuz measure /Z(m) is 
equal to the energy of M, i.e. 

e(M) = X -\ n {M) {dy). (15) 

o 

Therefore /i<M) is also called the energy measure of M. For M, L EM let 

(M,L) := I((M + L)-(M)-(L)). 

Then ((M, L) t )t>o is a CAF of bounded variation on each finite interval. Furthermore the 
finite signed measure [X{m,l) defined by H{m,l) '■= \{.H{m+l) — H{m) — /-*(£)) is related to 
(M, L) in the sense of (14). If G G B£, then J E Gd/i^.^ is symmetric, bilinear and positive 

o o 

on M x M- 
Define 

M c = {N\N is a finite CAF, e(N) = 0, £7 y [|JV t |] < oo 
for strictly S-q.e. y G E and all t > 0}. 

For F G JF, can uniquely be decomposed (see [11, Theorem 4.5. (i)], [13, Remark 0.17]) 

as 

A [F] = M [F] +N [F]^ M [F] eMj N lF] eA f c . (16) 
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The identity (16) means that both sides are equivalent as additive functionals w.r.t. 
Cap 1)giS . The uniqueness of (16) implies aM^ + bM^ = M^ aF+bG \ aN^ + bN^ = 
N [aF+bG]^ for any a ,beR, F,G eJ 7 . 
From Lemma 2.1 in [14] we know that for F G T 



a 2 



li( M [F])(dxds) = — (d x F) pdxds, 
and moreover, if F is constant pdy-a.e. on a Borel set B. Then 

V(MiF])(B) = 0. 

Now let us come back to the identification of Y. In order to identify the drift part we 
might proceed as follows. Denote by S x the Dirac measure in x G K. If j3(s) > —7(5) a.e. 
s, then integrating by parts in (7) we obtain that the generator of the diffusion is given 
informally in the sense of distributions by 

u 2 u 2 5 1 

LF(s, x) = —d xx F(s, x) + — - — ■ — —d x F(s, x) + d t F(s, x) + u F (dx)ds 

O O [X + 7(5)) 

where the boundary term vp is given by 

u F (dx)ds = [pd x F(s, x) — (1 — p)d x F(s, x)} — \x + 7(s)| <5_1 e ^~5^ s )(dx)ds 
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— I.T 4- - : > i: ' ' 



;i -p)<9+F(s,x)^|x + 7(s)| <5 ~ 1 e-^^ 7(s) (rfx)rfs, (17) 



and where d x , resp. 9^ , denote the right hand, resp. the left hand derivative in space. 

bx 2 

If we can show that \x + 7(s)| <5_1 e ~5 Ki ^(dx)ds, k : 1R + — > IR + locally bounded and 
measurable, is a smooth measure w.r.t. A, then there is a unique PCAF representing this 
measure by Theorem 2.2 in [14]. Theorem 2.3 in [14] then allows to identify the drift part. 
We will identify the corresponding diffusion when 5 > 1. 

Let M + x K. = lJn>i K-ni where (K n ) n >i be an increasing sequence of compact sub- 
sets of R + x ffiL Let E n :— K n C\ E, n > 1. Since Cq(E) C J 7 dense, it follows from 
[10, III. Remark 2.11] that (E n ) n >i is an £^-nest in the sense of [10, III. Definition 2.3(i)]. 
Consequently, P y (lim n _^ 00 a-^ < 00) = for £-q.e. y G E, hence in particular for pdy-a.e. 

y G E (see [10, IV. Lemma 3.10]). We obtain that (£' n ) n >i is an strict £-nest by (11). 
[13, Lemma 0.8 (ii)] now implies P 2/ (lim n _ >0O (Tgc < 00) = for strictly £-q.e. y G E. We 

may without loss of generality assume that E n C [0, n\ x ffi. n E, n > 1, and that E n is 
contained in the interior of E n+ \ for any n > 1. From now on we will fix such a strict 
£-nest (£ n ) n >i. 

A proof following similar lines of arguments to the proof of the next proposition can be 
found in [14]. However, we include it for the readers convenience. 
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Proposition 2.3 Let 5 > 1, k : IR + — > R &e measurable and locally bounded, such that 
(s, k(s)) G i? s /or eac/i s > 0. TTie measure 



Ie n (s,x)\x + ^(s^-'e-^S^^ds, N> 1, 
smooth w.r.t. (A,V). 

Proof We only show the statement for 5 > 1. The proof for 5=1 works in the same 
manner and is even easier since the derivative of y i— > 1 2/ 1 45 1 disappears as it is constant, 
so there are less additional terms (cf. below). Let N > 1. Let F G Cq(E), ip G C£°(R), 
< ^ < 1, l^-^loo < 2. Since « is locally bounded, the following two values 

KNmin ■= inf{/c(f)|t G [0,iV]}, /t7v max := sup{/«(t)|t G [0, AT]}, 

are finite. Let ip — 1 on [/cjv m j„, Kjvmax]) ^ = on [/tjvmax + 1, oof. For s G [0, AT] we have 



F(s,k(s))|k(s) + 7 (s)| 5 - 1 e — ^" = - / 9 X (^(a;)F(s,x)|x + 7(s)| 5 - 1 e-— ) cfe 



and thus 



^ E N 



\s-i. 1 



\F\(s,x)\x + >y(s)\ d 1 e 2 5 K{s) (dx)ds 



< 



N 



+1 



d x U(x)F(s,x)\x + 7(s)|' 5 " 1 e"^ N ) rfx 
< 2/ / + |F|)|a; + 7(s)| 5 - 1 e-^^s 

J k(s) 

N pK Nmax +l 



ds 



+ / / |F| ((S-l^x + ^s^-bxlx + ^s^-^e-^dxds 

J0 J k(s) 

< C N ^A £ 2(F,F) + I(F), (18) 
with 

1(F) := / / |F| ((5-l)|x + 7(s)| 5 - 2 -6a;|a; + 7(s)|^ 1 )e-^^s, 

and = f^r ( ; m " +1 |x + 7 ( S )| a - 1 e-^dxd S . 
Let AT C E be compact, and Cap" 4 (AT) = 0. By (8) 

Cap- 4 (AO = inf{A(F, F);Fe C^E)}, 

where Cl )K {E) = {F G C^(E)|F(s, x) > l,V(s,x) G AT}. Hence, there exists (F n ) neN C 
Cq(E), F n (s,x) > 1, for every n G N, (s,x) G AT, such that \F n \ v — > as ra — > 00. Since 
normal contractions operate on V we may assume that sup ngN sup^ ^^ |F n (s, x)| < C. 
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Selecting a subsequence if necessary we may also assume that Hindoo \F n \ = p(s, x)dxds- 
a.e, hence dxds-a.e. Consequently, using Lebesgue's theorem we obtain 

I(F n ) — > as n — > oo. 

Therefore by (18) 

f bx 2 

I l K (s,x)\x + ~i(s)\ 5 ' 1 5 K{s) (dx)ds 

J En 

< limsup _ \F n \(s,x)\x + *y(s)\ 5 ~ 1 e~ ! ^~5 K ( s )(dx)ds 

n— >oo J E N 

< limsup |c W y/^(F n ,F n ) + = 

Since 1^ (s, x)\x + 7(s)| <5 ~ 1 e i~5 K ( s )(da;)ds, as well as Cap" 4 are inner regular we obtain 



2 

that the measure (s, + 7(5) | 5_1 e i~ 5 K ( s )(dx)ds is smooth w.r.t. (^4, V). 

□ 



Let us choose {Jm)m>i, {Hm)m>\ C C 2 (£), with 

H M (s,x) : = 

M > 1, and 

J M {s,x) : = 

M > 1. Let further 
and 



x for (s, x) G -Em 
for (s,x)6^ +1 , 



s for (s, re) G -Em 
for (s, x) G E? M+l , 



H(s, x) := x, 
J(s, a;) := s. 



(-£/m)m>i (resp. (Jm)m>i) is a localizing sequence for /J (resp. J). Obviously 

4?4 = for an Y K > L > M. 

We claim that 

M™ c = c for any K > L > M. 

fc M fc M 

Indeed, for strictly £-q.e t/6£, and any t > 0, 
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E„ 



(M [Hk ~ Hl] ) 



tAcr-F 



= E„ 



< E„ 



tAo-^c 



fH M {Ys)d{M^- H ^) t 
Jo 



By Lemma 2.1(H) in [14] ^ i- M (y s )d{M\ H K-H L ] )s = h {m[ h k -h l]) (E m ) = 0. Thus by injectiv- 
ity of the Revuz-correspondence (see [11, Remark 5.2(h)]) E y jjj 1^ M (F S ) d{M^ Hl< - H ^)^ ■ 



strictly S-q.e. y G E. Hence the same is true for E y {M^ Hk Hl 



It Aa-F 



. We know that 



(M, 



[H K -H L U 2 



(Ml H * 



t>0 



is a martingale w.r.t. P y for strictly S-q.e. y G E. The 



optional sampling theorem then implies 



E„ 



[H K -H L U 2 



En 



(M [Hk ' Hl] ) 



tAa-E 



= 



for strictly E-q.e. y G E and the claim is shown. The analogous statements hold for A^ Jk \ 
M^ Jk \ Thus we may set 



M\ H] := lim M t 

M— >oo 



[%] 



and 

in order to obtain 



M\ A : lim m\ Jm] 

M^oo 



Nj* :=Af ] -M\ H \ 
N [ t J] := A[ J] - M [ /\ 



Af ] =Ml H] +Nl H \ 
[J] 



A [ t J] = M\ J] + N t 



[J] 



Note that n\ H] = lim M ^oo n\ Hm] , N l t J] = lim M ^oo M [ / Mi . We want to find the explicit 
expressions for M^ J \ N^ J \ M^ H \ N^ H \ Let F G Cq(E). Integrating by parts we obtain for 
any G G Cq{E) 



-£(F,G) 



^ ( 



— - — ^— - — bx J d x F J + d t F J Gpdxds 



2 roo /-/J(sj 2 

(1-p)- / / dJ^FG|z + 7 (s)| 5 -V^W S 
8 Jo y/3( s ) v y 



(7 



(^b)"^)^) + ^) G ^ 



+ - 



^ G9a.F{(l - p) I{^( s)> _ 7 ( s )} +pl{ /3 ( s )=- 7 ( s )}} |a; + 7(s)| <5 x e 6 2 <L 7 ( s) (cte)ds 



a 



+Y J E Gd x F(2p- l)I {/3(s)> _ 7W} |a; + 7(s)r~ 1 e--%)(rfa;)^. (19) 
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Obviously, (19) extends to G G V&. By Proposition 2.3, the measure 

2 

% M (s,a;)y|a; + 7(s)| 5 - 1 e-^(5 K(s) (rfa;)rfs, M > 1, 5 > 1, 

k = — 7,/9, is smooth w.r.t. (.4., V). Let £^ denote the unique positive continuous additive 
functional (PCAF) of F associated to ^-\x + / y(s)\ 5 ~ 1 e~^~5 K ( s )(dx)ds (see Theorem 2.2 

in [14]). Then f* G(Y s )d£g is associated to G(s, x)^-\x + 'y(s)\ 5 ~ 1 e~^~ 5 K ( s )(dx)ds for any 
G G B b (E). In particular, £ < ~ 7 vanishes, if 8 ^ 1. We obtain 

+(2p-l) / d x FI {0 

oJ>-7oJ} 

Jo 

+I {5=1} / 9 X F{(1 -p)I {/ 3oj>- 7 oj } +pI {/ 3oj=- 7 oj}} (Y s )d£r. (20) 

Indeed, if we denote the r.h.s. of (20) by A t then in particular by (19) — S(F, G{W) = 
lim^oo a 2 E diWpdy [/ °° e~ at A t dt] for all W G W 6 . Hence A"J F] = A t by Theorem 2.3 in 
[14]. On the other hand, by Lemma 2.1 (i) in [14] the Revuz measure /j>/mI f ]) is equal to 
^-(d x F) 2 pdy. A simple calculation shows that the Revuz measure of ^ J* * (d x F) 2 (Y s )ds is 
also equal to ^(d x F) 2 pdy. Consequently, we have (M^) t = ^ J* * (<9 :c F) 2 (F s )(is (see [11, 
Remark 5.2(h)]) and therefore we may assume that 

M\ F] = ^£d x F(Y s )dB s (21) 

with ((B t )t>o, P y , {J~t)t>o) being a Brownian motion starting at zero for strictly £-q.e. 
yEE. 

By (19), (20), applied to Jm, letting M — > 00, we obtain 

J(Y t ) = J(Y )+t. 

We put 

X t :=#(F t ), t>0, 

so that 

F t = (j(F ) + t,x t ). 

Applying again (19), (20), but this time to H M) letting M — > 00, we obtain 
X t = X + °-B t + £ f ~ bX s ds + (2p - 1) f I {/3oJ> _ 70j} (F s )^f 

+I {5 =i } / { (1 - p)I {/3 oj>- 7 oj} + pI {/3 oj=- 7 oj} } (Y s )d£;\ (22) 
Jo 
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(22) holds P r a.s for £-q.e. y E E. 



2.1.3 Identification of the process Y associated to the Generalized Dirichlet 
form involving the semimartingale local time 

The positive increasing processes £ -7 are so-called Markovian local times because they 
are PCAF's of a strong Markov process (namely Y) with related 1-potentials and associ- 
ated smooth measures that are singular (cf. (17), Proposition 2.3, and the Revuz formula 
(14)). It is a priori not clear whether these Markovian local times are semimartingale 
local times, i.e. increasing processes that occur in Tanaka's formula. There are Markov 
processes that are semimartigales and where the Markovian local times and the semi- 
martingale local times are different to each other (see e.g. [1]) but this will not be the 
case in our situation (see Remark 2.4(h)). 



Remark 2.4 (i) If 5 > 1, then (x + j(s)) 1 p(s,x)dxds is a smooth measure, since 
(x+'j(s))~ 1 E Lj oc (E , pdy) , and therefore X is a semimartingale. If 5 = 1, J* * — + ^)(y )) ^ s 
disappears, and X is again a semimartingale. 

If we had admitted 5 < 1, then clearly (x + , y(s))~ 1 ^ Lj oc (E,pdy), hence X would not be a 
semimartingale, and we would have to work with principal values in (22). Further Propo- 
sition 2.3 wouldn't apply. Nonetheless, it is clearly possible to consider the case 5 < 1 (see 
Remark 2.9 below). 

Finally, since all subsequent transformations applied to X (5 > 1!) are keeping the class of 
semimartingales invariant, the processes Y, Z, constructed below, and renamed as R, \[R, 
in the introduction, will remain semimartingales. 

(ii) In our case it will turn out that if k is regular enough, e.g. k E H^(M. + ), then 11 
(restricted to its support) is a constant multiple of a classical semimartingale local time 
(see e.g. [9, chapter VI] for a rigorous definition of semimartingale local time). We will 
determine these constants for k — f3, k — —7, but we remark that everything would have 
worked exactly in the same way up to (22) if we only assumed f3, —7, to be measurable and 
decreasing. The later Girsanov transformation of section 2.2 for which only 7 E H^(M. + ) 
is needed, can also in this case be applied to (22) (see Remark 2.8(H)). 



In order to determine the constants mentioned in Remark 2.4(h), we will evaluate the 
l.h.s. of Tanaka's formula (23) below with the help of (16) and then compare the result 
with the r.h.s. of (23) (see (25), (26)). We consider the point-symmetric derivative 




1 if x > 
if x = 
-1 if x < 0, 
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the left continuous derivative 

~sgn(x) 

and the right continuous derivative 

sgn(x) 



1 if x > 
-1 if x < 0, 



1 if x > 
1 if x < 0, 



of \x\. Let k be continuous, and locally of bounded variation. Since X — k is a continuous 
semimartingale, we may apply Tanaka's formula 

\X t - K (t)\ = \X -k(0)\+ [ f(X s -K(s))dX s + £° t f (X -«), (23) 

(cf. e.g. [9, VI.(1.2) Theorem, (1.25) Exercise]), where 

( £° t (X-K) if / = s^n 
£° t f (X-K) = \ f t + (X-K) if f = sgn 
( ^-{X-k) if f = sgn, 

and ^°(X — k) (resp. ^° + (X — /c), £®~{X — k)), is called the symmetric local time (resp. 
upper local time, lower local time) in zero of the continuous semimartingale X — k. In the 
book [9] they authors decided to work with the left continuous derivative of x i— > and 
they use the expression L\ for our £® + . For our framework it is more intuitive to work 
with the point symmetric derivative (see Remark 2.7 below). 
By [9, VI. (1.9) Corollary, (1.25) Exercise] we have P y -&.s. for S-q.e. y G E 

£°+(X - «) = Urn- / I m (X s - K(s))d(X,X) s , 
£° t -(X - K ) = lim- f\ iefi] (X a - K(s))d(X,X) s , 

and 

From [15] we know that {s > 0\X S — k(s) = 0} is P y -a.s. of Lebesgue measure zero for 
£-q.e. y e E. 



Let (F m )m>\ C Cq(E), such that 

F M (s,x) 



1 for (s, x) G -Em 
for (s,x)eF M+1 , 
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M > 1, and let 

— („ „\ — U _ „. ^ t 

L loc 



k(s,x) :— \x — k(s)\, k G if; 1. 



It is easy to see that kFm G J? 7 for any M > 1. We will use the same localization 
procedure as before. Thus, if Aff 1 := liniM^ A/f Fm] , and ivf ] := liniM^ TVf Fm] , then 
4 R 1 = Mf 1 + Tvf ] . If G G Cq 1 (£), then 

-£((3F M ,G) = — J f3d xx F M Gpdxds + — y 2sgn(j3)d x F M Gpdxds 

+ Y J e - + sgn@)F M ) Gpdxds 

-(1-p) y (Pd x F M - F M )GI Ws)> ^ l{s)} ^\x + -f(s)\ 5 ' 1 e- k ^5 (3{s) (dx)ds 

f a 2 b 2 

+p J {f3d x F M + F M )G—\x + 1 {s)\ s - 1 e- b -^5 f3(s) {dx)ds 

f - - a 2 bx 2 

-(1 -p) ((3d x F M - F M )sgn{(3)GI Ws ) > „ 1 ( s ) } — \x + -f(s)\ s ~ 1 e~^~8-^ s )(dx)ds 

+ / (\x - (3{s)\d t F M - F M sgn(x - (3(s)) (3\s))Gpdy. 
Je 

Obviously, the last equation extends to G G Thus, letting M — > oo, 



^f 1 = ^ 







if + 70 J 



- fctfj s<7n(/3)(y a )ds - jf S ^n(/3)(y s )d/3(J(y s )) 



+(1 - p) f I {poJ> _ loJ} (Y s )d£P + pf t 
Jo 

+{l-p) f sgn(P)(Ys)I{f]oj>^oj } (Y s )d£^. (24) 
Jo 

On the other hand by Lemma 2.1 (i) in [14] 

a 2 2 
/V^m!) = T dx ^ X ~ P( s >\ Fm > pdy 



a 2 



= -j ((\x- P(s)\d x F M y + 2\x- P(s)\d x F M sgn (x - (3(s)) F M ) pdy 



a 2 



+— (F M sgn(x - (3(s))) pdy. 



We obtain (AfW) t = \ f* sgn((3)(Y s ) 2 ds. Consequently, we may assume that 

Mf = ^j\gn(P){Y s )dB s . 
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Note that J* * sgn((3)(Y s )d£^ = 0, because for its associated signed smooth measure we 
have sgn(x — P(s))5 l 3( s - ) (dx)ds = 0, since sgn(0) = 0. Therefore 

\X t -(3(J(Y t ))\ = \X -(3(J(Y ))\+ f sgn{X s -(3{J{Y s )))d{X s -(3{J{Y s ))) 

Jo 

+Kf + (1 - P) I I { ^oj>- 7 oj}(F s )^. (25) 
Jo 

For and?iGR+ define 

K u (t) := k{u + i). 

Recall that 

P y (J(Y t ) = J(y) + t) = l, 

so that 

K J(y)i t ) = K ( J ( Y t)) Py-a-s. 

Comparing (25) with Tanaka's formula (23) we see that 

f t (X - (3 J{y) ) = pi? + (1 - p) f I {/3j(y)(s)> _ Tj(j , )(s)} ^ P.-a.s (26) 

Jo 

for S-q.e. y G E. In a similar way one can see that we have 

£° t (X + u {y) ) = \p t + {X + lJ{y) ) 

= / P I {A/( y) ( S )=-7.7( !/ )W} + i 1 -p) I {A/ (y) ( S )>-7./( !/ )W} rf C 7 , (27) 
Jo 

Py-a.s. for S-q.e. y G E. Therefore, (22) rewrites P^-a.s. as 

Xt = X Q + a -B t + ( ^ f 5 ~ l -bX s ds 
2 8 i X s + 7j ( j /) (s) 

+ (2p - 1) jf Itf J(w)W >- 7j{s) «}<(* - /3 J( , } ) + + 7JW ). (28) 

2.2 Girsanov transformation of (4) and concluding remarks 

Let b G R + . We define 

1 /"* 

J 
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and 

d n = e -^/o^j( B )W+^(«)W^--3^^l^(»)W+^^(»)Wl a<is dp on t, 

v y y t' 

Obviously, Novikov's condition is satisfied since 7 G /// 'f(M + ), hence is a Brownian 
motion under the equivalent measure Q y . Put 

Y t :=X t + lj(Yo) (t). 

Then 

- (3 J{y) ) =i° t (X + 1Jiy) - {(3 J{y) + lJ{y) )) =i° t (Y- X J(y) ) Q,-a.s, 

since Q y is equivalent to P y . Analogously, ^ + (X + 7j( y )) = ^ + {Y) holds Q y -a.s. Thus, 
under Q y , Y t > (since X t > — 7(£)), and 

+ (2p - 1) " W + %^? + (n (29) 

where QyC^o = + l(J(y))) — 1- Moreover, under Q^, Z t := F t 2 , satisfies 

^ = ^0 + o- / v^dW, + I" / ( 5 " 6Z ^ S 
Jo 4 </o 

+(2p-l) [ I {Aj(y)(s)>0} 2 v / ^<(v / Z-A J( , ) ), (30) 
and Z = (#(?/) + 7 (J(y))) 2 Q„-a.s. 

Remark 2.5 Setting 6 = 0, and replacing a 2 b r yj( y )(s) by a 2 c, c G m £/ie expression 
for Q y , we construct instead of (30) a positive solution Z = Y 2 to 

Z t = Z + a [ ^Z~ s dW s + ^ / (5- c^Z~ s )ds 
Jo 4 Jo 

+(2p-l) [ I {Xj(y)is)>0} 2^ s d£ s (^Z-X Jiy) ). (31) 
Jo 

For p = |, zs well-known as the double square-root (DSR) model of Longstaff in 

financial mathematics. Forp = | ; (30) is the well-known Cox-Ingersoll-Ross model (CIR). 

Remark 2.6 The equations (28), (30), (31), are in the sense of equivalence of additive 
functionals. This means that they hold for initial conditions outside some exceptional set. 
If X 2 is constant, say X 2 = c > 0, we are in the symmetric case, and it is clear that 
the parabolic capacity is comparable with the elliptic one, so that we can start from every 
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X , Y , Z = x > 0, if 5 < 2, and for every X , Y , Z = x > 0, if 5 > 2. Indeed Y is then 
associated to the symmetric (!) Dirichlet form (cf. [3]), which is uniquely determined as 
the closure of 

poo 2 

£ p (f,g):=J o -xf\x)g\x)xi- l e-^p(x)dx- f,g& C^°([0, oo)) 

in L 2 ([0, oo), x^e~^ p(x)dx), where p(x) = (1 — p)I{ x < c } + pl{ x > c j, and it is clear that 
the capacities are all equivalent for p G (0,1). p = \ corresponds to case without skew 
reflection. 

The regularity of the equations (28), (30), (31), i.e. the question whether we can start 
pointwise in the parabolic case, i.e. if X 2 ^ const is subject of forthcoming work. Note 
however, that the structure of (28), (30), (31), is not influenced by these questions of 
regularity. 

Remark 2.7 (i) Rewriting (25) with ~sgn,sgn, one can easily see that ^ + (X — j3) — 
pit , and i° t ~(X — (3) = (1 — p) J* I{ / g( s ) > _ 7 ( s )}d£f. This implies the following relations for 
/'/(/.'• A 2 ) in II}: 

£° t (R - A 2 ) = - A 2 ) = ^pjtt'iR - A 2 )- (32) 

Note that a stronger version of (32) holding for more general A 2 is obtained in [15] by 
purely probabilistic methods. One observes immediately the discontinuity of the local times 
in the space variable, thus we provide another example of diffusion with discontinuous 
local time (see e.g. [17]). Moreover, if \p\ > 1, then any of these local times is identically 
zero. Consequently, the associated process is the CIR process, for which uniqueness in any 
sense is known to hold. Consider the time dependent Dirichlet form £ (cf e.g. [7], [14]) 
corresponding to the time dependent CIR process (t,R t ), then 

poc poo g.2 

S(F,G) := / / —\x\d x F(t J x)d x G(t,x)\x\^~ 1 e'^dxdt 
Jo Jo 2 

poo poo 

II 5 i bx 

- / d t F(t,x)G(t,x)\x\^- 1 e~^dxdt. 
Jo Jo 

One can easily see that the local time £°(R — A 2 ) is uniquely associated to the measure 

— \x\2e~^ 5 X 2(t)(dx)dt. 

But this measure doesn't vanish if A 2 is different from zero on a set of positive Lebesgue 
measure. Therefore £°(R — A 2 ) cannot vanish identically, a contradiction. Thus an asso- 
ciated process cannot exist if \p\ > 1. 

(ii) The first reason to employ point symmetric derivatives, is that these better work out 
the intuitive structure of skew reflection, namely 2p — 1 = p — (1 — p), so upper reflec- 
tion with probability p, and lower reflection with probability 1 — p. At least in the case 
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(5 = 1, o" = 2,6 = 0, A 2 = 0) the just mentioned intuitive structure is rigorously described 
for the square root process (the skew BM) through excursion theory (see e.g. [6], [17]). The 
second reason is that symmetric local times correspond to symmetric derivatives, which 
are used in distribution theory, and therefore correspond to our analytic construction of 
the Markov process generator. 



Remark 2.8 (i) The relation (3) can easily be derived by writing down Fukushima's 
extended decomposition (16) in localized form for \Z t — \ 2 (t)\ = \(X + lit)) 2 — A 2 (t)| (Z 
as in (30)) and then comparing it with the symmetric Tanaka formula. More precisely, 
we obtain 

2^Z~ s dl s (VZ-\ J(y) ) = d£° s (Z - \ 2 J{y) ) 

by doing so. 

(ii) Coming back to Remark 2.4(H), suppose that we had assumed A = f3 + 7, where 
only 7 G H^ 2 (R + ), and (3 not necessarily continuous, but decreasing. Then, we would 
have obtained (29), (30), except that £°{\fz — Xj( y )) has to be replaced by £ x , where £ x is a 
positive continuous additive functional ofY, which only grows when Y = \, or equivalently 
Z = \ 2 . 

Remark 2.9 (The case 5 G (0, 1)) Since we are no longer in the semimartingale case 
we can no longer make use of the Girsanov formula with 7. Therefore one puts 7 = 0. 
Then one may still start as before with (7) and the following density 

bx 2 

p(t,x) = ((1 - p)I [0Mt)) (x) +pI[A(t),oo)(aO) M^V^. 

Note that p(-,x) is still assumed to be increasing in t, and that Proposition 2.3 could no 
longer be available. 

One may preferably directly start with the squared process since the technique of changing 
the measure will not be used. Thus one could proceed with the following (cf. (10)) time 
dependent form 

poo poo 2 

S(F,G) = —\x\d x F(s,x)d x G(s,x)p(s,x)dxds 
Jo Jo 2 

d t F(s, x)G(s, x)p(s, x)dxds, 



with 

p(t,x) = ((1 -p)I[o,A a (t))(aO +Z% 2 (t),oo)(aO) |x|2 _1 e~^\ 

and where p G (0,1), and A 2 G H^l(M, + ), are chosen, such that p is increasing in t. 
In order to convince the reader, we just line out the argument for the existence of the 
corresponding local time on A 2 . In fact one only has to show that 

— \x\*e~~*8xnt){dx)dt 

is smooth with respect to the symmetric part of £ . This can be done analogously to Propo- 
sition 2. 3. Of course, if A 2 is a constant, we consider the symmetric Dirichlet form of 
Remark 2. 6. 
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